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A quantum computer is a device that processes in- 
formation in a quantum-mechanically coherent fash- 
ion Era. In principle, it can exploit coherent quantum 
interference and entanglement to perform computa- 
tions, such as factoring large numbers or searching an 
unsorted—database, more rapidly than classical com- 
puters Elnuu. Noise, decoherence, and manufactur- 
ing problems make-constructing large-scale quantum 
computers difficult m 13 l. Ion traps and optical, cavities 
offer promising experimental approaches Lj'Ej, but no 
quantum algorithm has yet been implemented with 
those systems. On the other hand, because of their 
natural isolation from the environment, nuclear spins 
are particularly good 'quantum bits' E£l and their use 
for quantum computation is possible by applying nu- 
clear magnetic resonance. (NMR) techniques in an un- 
conventional manner I 1 Tl 19 l. Here, we report on the 
experimental realization of a quantum algorithm us- 
ing NMR, to solve a purely mathematical problem 
in fewer steps than is possible classically. In partic- 
ular, our simple quantum computer can determine 
global properties of an unknown function using fewer 
function 'calls' than is possible using a classical com- 
puter. 

We implemented the simplest possible version of the 
Deutsch-Jozsa (D-J) quantum algorithm u, which determines 
whether an unknown function is constant or balanced. A con- 
stant function f(x) from N bits to one bit either has output 
f(x) = for all x, or f(x) = 1 for all x. A balanced function 
has f(x) = for exactly half of its inputs, and f(x) = 1 for 
the remaining inputs. To determine with certainty whether a 
function is constant or balanced on a deterministic classical 
computer, requires up to 2 JV_1 + 1 function calls: even if one 
has looked at half of the inputs and found f(x) = for each, 
one still can't conclude with certainty that the function is 
constant. Injajintrast, the D-J algorithm, as improved by R. 
Cleve, et al. B3 and Alain Tapp, allows a quantum computer 
to determine whether f(x) is constant or balanced using only 
one function call. 

The D-J algorithm is well illustrated by its simplest possible 
case, when / is a function from one bit to one bit; this is the 
version that we have realized (it is also the simplest instance 
of Simon's algorithm 13). There are four possible f's, two of 
which are constant, fi(x) — 0, fz(x) = 1 and two of which 
have an equal number of and 1 outputs: feix) = x, f4(x) = 
NOT x. To determine whether such a function is constant or 
balanced is analogous to determining whether a coin is fair - 
with heads on one side and tails on the other; or fake - with 
heads on both sides. Classically, one must look at the coin 
twice, first one side then the other, to determine if it is fair 



or fake. The D-J algorithm exploits quantum coherence to 
determine if a quantum 'coin' is fair or fake while looking at 
it only once. The algorithm requires one 'input' spin and one 
'work' spin, and is schematically represented by the quantum 
circuit shown in Fig. |l| 

Experimentally, this quantum algorithm was implemented 
using the nuclear spins of the 1 H and 13 C atoms in a carbon- 13 
labeled chloroform molecule (CHCI3) as the input and work 
quantum bits ('qubits'). |0) (|1)) describes the spin state 
aligned with (against) an externally applied, strong static 
magnetic field Bo in the +z direction. The reduced Hamilto- 
nian for this 2-spia system is to an excellent approximation 
given by (h = 1) PA 



H 



-ujaIzA — ubIzb + 2-kJIzaIzb + T~i e 



The first two terms describe the free precession of spin A ( 1 H) 
and B ( 13 C) about -B with frequencies uj a /2tt « 500 MHz 
and lob/^tv ~ 125 MHz. I z a is the angular momentum op- 
erator in the +£ direction for A. The third term describes 
a scalar spin-spin coupling of the two spins of J « 215 Hz. 
~Henv represents couplings to the environment, including inter- 
actions with the chlorine nuclei, and also higher order terms 
in the spin-spin coupling, which can be disregarded (as will 
be described below). 
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FIG. 1. Quantum circuit for performing the D-J algo- 
rithm. (TO) Start with both the 'input' and 'work' qubits 
(A and B) in the state |0). (Tl) Perform the transforma- 
tion F:|0) - (|0> + |1»/V2, |l>_-> ( - |0> + |1»/V2, to A, 
and the inverse transformation Y to B, resulting in the state 

\ ELo (1^) — 1 1))" ^ ne input qubit in some quantum 
sense registers both and 1 at once. (T2) Call the function: 
apply / to A, and add the result to B modulo 2. As long as 
the quantum logic operations needed to evaluate / are carried 
out coherently, the work qubit now contains in some quan- 
tum sense the outputs of / on all possible-jnputs, an effect 
that Deutsch termed 'quantum parallelism' a. The two qubits 
are now in the state ^^2l, =0 \x) (\0 + f(x)) — \1 + f(x))) 

= sELoC -1 )^'! 31 ) (I ) - I 1 ))- ( T3 ) Perform the inverse 
of the transformations of (Tl), thereby taking the qubits out 
of their superposition states. If / is constant, then the fac- 
tors (—l)f' x ' are either all +1 or all —1, and the result of the 
transformation in this step is the state ±|00). If / is balanced, 
then exactly half of the factors ( — l)-^w are +1 and half are 
— 1, and the result of the transformation is the state ±|10). 
(T4) Read out A. If it is 0, then / is constant; if it is 1, then 
/ is balanced. 

The five theoretical steps of the quantum algorithm, (TO)- 
(Tl), were experimentally implemented as follows: 
(E0) An input state is prepared with a 200 mM, 0.5 ml sam- 
ple of chloroform dissolved in d6-acetone, at room tempera- 
ture and standard pressure. The O(10 18 ) molecules in this 
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bulk sample can be thought of as being independent single 
quantum computers, all functioning simultaneously. The the- 
oretically ideal result is obtained when the spins in all the 
molecules start out in the 00 state. Because the experiment is 
performed at room temperature, however, the initial density 
matrix p for the thermally equilibrated system has popula- 
tions diag(p) = [noo, noi, nio, mi] in the 00, 01, 10, and 11 
states, respectively, where p is the density matrix, and rii are 
proportional to e' Ei/kT /2 N w (1 - E t /kT)/2 N , with E t the 
energy of state i and N — 2 being the number of qubits used in 
our experiment. A variety of techniques exist to extcacLJrom 
this thermal state just the signal from the Il£j state [U'lj; we 
adopted the method of 'temporal averaging' E3, which involves 
the summation of three experiments in which the populations 
of the 01, 10, and 11 states are cyclically permuted before per- 
forming the computation. The essential observation is that 
[noo, noi, nio, nil] + [noo,nn,noi,nio] + [noo,nio,nn,noi] = 
a[l, 1,1,1]+ S[l, 0, 0, 0], where a — noi + nio + nn is a back- 
ground signal which is not detected, and 5 = 3noo — a is a 
deviation from the uniform background whose signal behaves 
effectively like the desired pure quantum state, 1 00) . The per- 
mutations are performed using methods similar to those used 
for the computation, described next. This technique avoids 
the technical difficulties of detecting the signal from a single 
nuclear spin, and allows a sample at room temperature, which 
produces an easily detectable signal, to be used for quantum 
computation. 

Note that while this method requires f(x) to be evaluated 
3 times, it is actually not necessary. Although step (TO) stip- 
ulates an input pure state 1 00) , the algorithm works equally 
well if the input qubit is initially |1); furthermore, when the 
work qubit is initially jl), it fails, and cannot distinguish con- 
stant from balanced functions, but this does not interfere with 
other computers which have worked. Thus, a thermal state 
is a good input for this algorithm, and only one experiment 
needs to be performed. Data from both thermal and pure 
state inputs are presented below. 

(El) Pulsed radio frequency (RF) electromagnetic fields are 
applied to transform the qubits as prescribed in (Tl). These 
fields, oriented in the x — y plane perpendicular to Bo, selec- 
tively address either A or B by oscillating at frequency uja or 
ujb- Classically, an RF pulse along y (for example) rotates a 
spin about that axis by an angle proportional to w tP, the 
product of the pulse duration t and pulse power P. In the 
'bar magnet' picture, a tv/2 pulse along y (we shall call this 
Y ) causes a z oriented spin to be rotated by 90° , onto x (sim- 
ilarly, we shall let Y denote n/2 rotations about —y, and X 
denote n/2 rotations about x, and so forth; subscripts will 
identify which spin the operation acts upon). This descrip- 
tion of the state is classical in the sense that a bar magnet 
always has a definite direction. In reality, however, a nuclear 
spin is a quantum object, and instead of being aligned along 
x, it is actually in a superposition of being up and down, 
(|0) + |1»/a/2. Likewise, a spin classically described as being 
along —x is actually in the state (|0) — \l})/\/2. (El) thus 
consists of applying the two RF pulses YaYb- 
(E2) The function y — > y © f(x) is implemented using RF 
pulses and spin-spin interaction. Recall that spin A represents 
the input qubit x, and B the work qubit y where / stores its 
output, /i is implemented as t/2 — XbXb —t/2 — XbXb, to 
be read from left to right, where t/2 represents a time interval 



of 1/4 J w 1.163 ms, during which coupled spin evolution 
occurs. Dashes are for readability only, and typical pulse, 
lengths were 10-15 ps. This is a well known refocusing E3 
pulse sequence which performs the identity operation, /a is 
t/2 — XbXb — t/2, similar to f 1 but without the final pulses, 
so that B is inverted, fz is Yb — T — YbXb — YaXaYa, 
which implements a 'controlled-NOT' operation, in which B is 
inverted if and only if A is in the |1) state. The naive 'bar 
magnet' picture can be used to get a feeling for how this works 
in case the inputs are 00 or 10, for which the subsequence 
Yb — t — Xb suffices (note that after (El), both spins are not 
just |0) or |1) but in a superposition of both, in which case 
the extra pulses of fz are necessary ED). First, Yb rotates B 
to +x. B then precesses in the x — y plane, about —i. Due 
to the spin-spin coupling, B precesses slightly slower (faster) 
if A = (A = 1). After r seconds, B reaches +y (— y) in the 
rotating frame. Xb then rotates B to +z (—2), i.e. to or 
1, where the final state of B depends on the input A. The 
precise quantum description is easily obtained by multiplying 
out the unitary rotation matrices. Finally, is implemented 
as Yb—t—YbXb — YaXaYa, which is similar to fz but leaves 
B inverted. 

(E3) The inverse of (El) is done by applying the RF pulses 
YaYb to take both spins back to ±2. Spin A, which was |0) 
at the input, is thus transformed into |0) or jl) for constant 
or balanced functions respectively. 

(E4) The result is read out by applying a read-out pulse Xa 
to bring spin A back into the x — y plane. The time varying 
voltage V(t) induced by the precession of spin A about B 
is recorded by a phase sensitive pick-up coil. Inspection of 
the spectrum of V(t) after a single experiment run and an 
appropriate read-out pulse, immediately reveals whether f(x) 
is constant or balanced, as shown in Fig. ^. 

We also characterized the entire deviation density ma- 
trix pA = p — Tr(p)//4 (Fig. ^) describing the final 2-qubit 
state. These results unambiguously demonstrate the complete 
proper functioning of the quantum algorithm, and provide 
data for the error analysis described below. 

Quantum computation requires that a coherent superposi- 
tion be preserved for the duration of the computation. This 
requires a highly isolated quantum system (small 7i en v), and 
fortunately, nuclear spins are naturally well-isolated from 
their environment. Phase randomization due to Bo inhomo- 
geneities was minimized by using about 30 electromagnetic 
coils to shim the static field to be constant to about one part 
in 10 9 over the sample volume. The longitudinal and trans- 
verse relaxation time constants Ti and T2 were measured us- 
ing standard inversions-recovery and Carr-Purcell-Meiboom- 
Gill pulse sequences E3, giving T\ w 19 and 25 seconds, and 
T2 ~ 7 and 0.3 seconds, respectively, for proton and car- 
bon; these were much longer than required for our experiment, 
which finished in about 7 milliseconds. 

The single most important source of errors in the experi- 
ments was the RF field inhomogeneity and pulse length cali- 
bration imperfections. A direct measure of this inhomogene- 
ity is the « 200 ps time constant of the exponentially de- 
caying envelope observed from applying a single pulse, as a 
function of pulse width. Including the population permuta- 
tion sequence, about 7 pulses are applied to each nucleus, with 
a cumulative duration of « 70 — WOps. 
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FIG. 2. Proton spectrum after completion of the D-J 
algorithm and a single read-out pulse Xa, with an ef- 
fectively pure input state 1 00) and with a thermal input 
state [Inset]. The low (high) frequency lines correspond 
to the transitions |00) 1 10) (|01) <-> [11». The fre- 
quency is relative to 499755169 Hz, and the amplitude 
has arbitrary units. The spectrum is the Fourier trans- 
formed time varying voltage V(t), induced in the pick 
up coil by the precession of spin A about B , at fre- 
quency loa, after the read-out pulse Xa- V(t) is given by 



V(t) ~ V Tr[e- 



^p(0)e l 



■IvjHt x ( . 



i&xA 



(TyA) 



where &{ X:V } are Pauli matrices, and p(0) is the density ma- 
trix of the state immediately before the readout pulse. By 
this convention, a spectral line for spin A is real and posi- 
tive (negative) when spin A is |0) (|1)) right before the Xa 
read-out pulse. Experiments were performed at Stanford 
University using an 11.7 Tesla Oxford Instruments magnet 
and a Varian UNIT Inova spectrometer with a triple-resonance 
probe. 13 C-labeled CHC1 3 was obtained from Cambridge Iso- 
tope Laboratories, Inc. [CLM-262] . 

The second most important contribution to errors is the low 
carbon signal-to- noise ratio, signal peak height/RMS noise w 
35, versus w 4300 for proton. The carbon signal was much 
weaker because the carbon gyromagnetic ratio is 4 times 
smaller, and the carbon receiver coil is mounted more re- 
motely from the sample. Smaller contributions to errors came 
from incomplete relaxation between subsequent experiments, 
carrier frequency offsets, and numerical errors in the data 
analysis. 

For this small-scale quantum computer, imperfections were 
dominated by technology, rather than by fundamental issues. 
However, NMR quantum computers larger than about 10 
qubits will require creative new approaches, since the sig- 
nal strength decays exponentially with the munber of qubits 
in the machine, using current schemes oEEt for N spins 
the signal from the initial state 00 ... is proportional to 
H00...0 oc NZ~ N where the single spin partition function 
Z ~ 2 at high temperatures. Furthermore, coherence times 
typically decrease for larger molecules, while the average logic 
gate duration increases. Nevertheless, there is hope; for ex- 
ample, due to the ensemble nature of the NMR approach, one 
can infer the output result as long as a distinguishable major- 



ity of the molecules reaches the correct final state. Creating 
an effective pure state is thus not always necessary, as we 
have demonstrated. Optical pumping and other cooling tech- 
niques can also be used to pre-polarize the sample to increase 
the output signal amplitude, since Z m 1 at low temperatures. 
Quantum computation clearly poses an interesting and rele- 
vant experimental challenge for the future. 
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FIG. 3. Experimentally measured and theoretically ex- 
pected deviation density matrices after completion of the D-J 
algorithm. The diagonal elements represent the normalized 
populations of the states 1 00) , 1 01) , 1 10} and 1 1 1} (from left 
to right). The off-diagonal elements represent coherences be- 
tween different states. The magnitudes are shown with the 
sign of the real component; all imaginary components were 
small. The deviation density matrix was obtained from the 
integrals of the proton and carbon spectral lines, acquired 
for a series of 9 experiments with different, read-out pulses 
for each spin (quantum state tomography EJ). The observed 
experimental non-idealities can be quantified as follows. In 
the experiments, the normalized pure-state population (ide- 
ally equal to 1), varied from 0.998 to 1.019. The other 
deviation density matrix elements (ideally 0), were smaller 
than 0.075 in magnitude. The relative error e on the ex- 
perimental pure-state output density matrix p exp , defined as 
e =|| Pexp — ptheory \\ / \\ ptheory ||, varied between 8 and 12%. 

Note: during the preparation of this manuscript we became 
aware of a closely related experiment by J. A. Jones and M. 
Mosca at Oxford University EE 



1. Deutsch, D. 'Quantum Theory, the Church- Turing Principle and 
the Universal Quantum Computer,' Proc. R. Soc. Lond., A, 400, 
97-117 (1985). 

2. Shor, P., 'Algorithms for quantum computation: discrete loga- 
rithms and factoring,' Proc. 35 th Ann. Symp. on Found, of Com- 
puter Science (IEEE Comp. Soc. Press, Los Alomitos, CA, 1994) 
124-134 . 

3. Divincenzo, D. P., 'Quantum computation,' Science 270, 5234, 
255-261 (1995). 

4. Lloyd, S., 'Quantum-mechanical computers,' Scientific American 
273, 44-50 (Oct. 1995). 

5. Ekert, A. and Jozsa, R., 'Quantum computation and Shor's fac- 
toring algorithm,' Rev. of Mod. Physics 68, 3, 733-753 (1996). 

6. Deutsch, D. and Jozsa, R., 'Rapid solution of problems by quan- 
tum computation,' Proc. R. Soc. Lond.. A 439, 553-558 (1992). 

7. Simon, D., 'On the power of quantum computation,' Proc. 35 th 
Ann. Symp. on Found, of Computer Science (IEEE Comp. Soc. 
Press, Los Alamitos, CA, 1994) 116-123. 

8. Grovcr, L. K., 'Quantum computers can search arbitrarily large 
databases by a single query,' Phys. Rev. Lett. 79, 23, 4709-4012 
(1997). 

9. Unruh, W. G., 'Maintaining coherence in quantum computers,' 
Phys. Rev. A 51, 2, 992-997 (1995). 



3 



10. Chuang, I. L.. Laflammc, R., Shor, P. and Zurck, W. H., 'Quan- 
tum computers, factoring, and dccohcrcncc,' Science 270, 5242, 
1633-1635 (1995). 

11. Landaucr, R., 'Dissipation and Noise Immunity in Computation 
and Communication,' Nature, Vol. 335, 779-784 (1988). 

12. Landaucr, R., 'Is Quantum Mechanics Useful?' Phil. Trans. R. 
Soc. bond. A, 335, 367-376 (1995). 

13. Palma, CM., Suominen, K.-A. and Ekert, A.K., 'Quantum Com- 
puters and Dissipation,' Proc. R. Soc. Lond., A, 452, 567-584 
(1996). 

14. Monroe, C, Meekhof, D.M., King, B.E., Itano, W.M. and 
Wincland, D.J., 'Demonstration of a fundamental quantum logic 
gate,' Phys. Rev. Lett. 75, 25, 4714-4717 (1995). 

15. Turchette, Q. A., Hood, C.J., Langc, W., Mabuchi, H. and Kim- 
ble, H.J., 'Measurement of conditional phase shifts for quantum 
logic,' Phys. Rev. Lett. 75, 25, 4710-4713 (1995). 

16. Lloyd, S., 'A potentially realizable quantum computer,' Science 
261, 5128, 1569-1571 (1993). 

17. Gershenfeld, N. and Chuang, I. L., 'Bulk spin-resonance quantum 
computation,' Science 275, 5298, 350-356 (1997). 

18. Cory, D. G., Price, M. D., Fahmy A. F. and Havel, T. F., 'Nu- 
clear magnetic resonance spectroscopy: an experimentally accessi- 



ble para digm for quantu 



E-print |quant-ph/9709001 
19. Cory, D. G., Fahmy, A. 



computing,' Physica D, in print; LANL 
gov (1997). 

and Havel, T. F., 'Ensemble quan- 



tum computing by NMR spectroscopy,' Proc. Nat. Acad. Sci. 94 
1634-1639 (1997) 

20. Clcve, R., Ekcrt, A., Macchiavello, C. and Mosca, M.,, Proc 



R , Soc . Lo pd., A, 454, 339-354 (1998): LANL E-print juant 



ph/9708016 
21. Slichter, c' 



P., Principles of Magnetic Resonance. (Springer, 
Berlin, 1990). 

Knill, E., Chuang, I. L. and Lafiamme, R. , 'Effective pure states 
for bulk quantum co input ation , ' — U* — aapear in Phys. Rev. A, 
(1998); LANL E-print 



quant-ph/9706053 



23. Ernst, R. R., Bodenhauscn, and Wokaun, A., Principles of Nu- 
clear Magnetic Resonance in One and Two Dimensions. (Oxford 
University Press, Oxford, 1994). 

24. Chuang, I. L., Gershenfeld, N., Kubinec, M. G. and Leung, D. W., 
'Bulk quantum computation with nuclear magnetic resonance: 
Theory and experiment,' Proc. R.. Soc. Lond., A, 454, 447-467 
(1998). 

25. Warren, W. S., 'The Usefulness of NMR Quantum Computing,' 
Science, 277 1688-1690 (1997). 

26. Jones, T. F. and Mosca, M., 'Implementation of a Quantum Al- 
gorithm to Solve Deutsch's Problem on a Nuclear Magnetic Res- 



onance Quanti 
LANL E-print 



m Hnmpntnr ' cnl 



quant-ph/9801027 



to J. of Chem. Phys. (1998); 



Acknowledgments 



We thank Alex Pines and Mark Kubinec for helpful discussions. 
This work was supported by DARPA under the NMRQC initiative. 
L.V. gratefully acknowledges a Francqui Fellowship of the Belgian 
American Educational Foundation and a Yansouni Family Fellow- 
ship. 

Correspondence and requests for materials should be addressed 
to I.L. Chuang, electronic address ichuang@almaden.ibm.com 



4 



